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Abstract 

Riemann theta functions are used to construct one-periodic and two-periodic 
wave solutions to a class of (2 + l)-dimensional Hirota bilinear equations. The 
basis for the involved solution analysis is the Hirota bilinear formulation, and the 
particular dependence of the equations on independent variables guarantees the 
existence of one-periodic and two-periodic wave solutions involving an arbitrary 
purely imaginary Riemann matrix. The resulting theory is applied to two nonlin- 
ear equations possessing Hirota bilinear forms: ut + Uxxy — Suuy — 3uxV = and 
ut + Uxxxxy — {5uxxV + lOuxyU — 15u^v)x = where Vx = Uy, thereby yielding their 
one-periodic and two-periodic wave solutions describing one dimensional propagation 
of waves. 
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1 Introduction 



It is always important to search for exact solutions to nonlinear differential equations. 
Different approaches, particularly in soliton theory, provide many tools for constructing 
explicit and exact solutions. Various kinds of exact solutions such as solitons, positons, 
complexitons, solitonoffs and dromions have been presented for nonlinear integrable equa- 
tions Successful methods include the inverse scattering transform the Darboux 
transformation [2j, Hirota direct method [3], and algebro-geometrical approach [1]. 

The algebro-geometrical approach presents quasi-periodic or algebro-geometric solu- 
tions to many soliton equations, which contain the KdV equation, the sine-Gordon equa- 
tion and the nonlinear Schrodinger equation. In recent years, such an approach have been 
applied to many (2-|-l)-dimensional nonlinear integrable equations [in]-[I3]- Nonlineariza- 
tion of Lax pairs [I1]-[I7] plays a crucial role in connecting the resulting algebro-geometric 
solutions with Liouville integrable Hamiltonian systems. The approach, however, needs 
Lax pair representations and involves complicated calculus on Riemann surfaces. 

On the other hand, the Hirota direct method provides a powerful way to derive soliton 
solutions to nonlinear integrable equations and its basis is the Hirota bilinear formulation 
[3]. Once the corresponding bilinear forms are obtained, multi-soliton solutions and rational 
solutions to nonlinear differential equations can be computed in a quite systematic way 
[3], even through Wronskian, Casoratian or Pfaffian determinants [18]- [M]- It is based 
on Hirota bilinear forms that Nakamura presented an approach to multi-periodic wave 
solutions of nonlinear integrable equations [25], together with applications to various typical 
examples of integrable equations. Such a method of solution does not need any Lax pairs 
and their induced Riemann surfaces for the considered equations. The presented multi- 
periodic solutions can be reduced to soliton solutions under asymptotic limits [261 EZ] • The 
advantage of the method is that it only relies on the existence of Hirota bilinear forms. 
Moreover, all parameters appearing in Riemann matrices are completely arbitrary, whereas 
algebro-geometric solutions involve specific Riemann constants, which are usually difficult 
to compute. 

In this paper, motivated by Nakamura's idea [25], we would like to use Riemann theta 
functions to generate one-periodic and two-periodic wave solutions to a particular class of 
(2-1-1 )-dimensional Hirota bilinear equations, and the corresponding solution analysis will 
be made to guarantee the existence of one-periodic and two-periodic wave solutions to the 
selected class of (2-|-l)-dimensional nonlinear equations. As illustrative examples of the 
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resulting theory, we will discuss two nonlinear equations possessing Hirota bilinear forms: 

Ut + Uxxy — 'iuuy — ?)UxV = and Ut + Uxxxxy — i^UxxV + lOuxyU — 15u'^v)x = 0, 

where Vx = Uy, and their one-periodic and two-periodic wave solutions involving an arbi- 
trary purely imaginary Riemann matrix will be explicitly presented. 

2 Existence of one-periodic and two-periodic wave so- 
lutions 

Let us consider an evolution equation in 2 + 1 dimensions: 

Ut = K{u,Ux,Uy,- ■ (2.1) 

where t G M is the time variable and x,y ^ M. are the space variables. We assume that 
under a transformation 

u = Uo-2(lnf{x,y,t))xx, (2.2) 

where Uq is a special solution to 02. ip . the evolution equation (12.11) can be transformed into 
a Hirota bilinear equation 

F{Dx,Dy,Dt)f-f = 0, (2.3) 

where F is a polynomial in the three variables. Here and below, the Hirota bilinear 
differential operators [3J are defined by 

DPD'^D'J{x,y,t)-g{x,y,t) 

(2.4) 

= {dx - dx'Yidy - dy>Y{dt - dt'Yf{x, y, t)g{x', y', t')\x'=x,y'=y,t'=t, 

where p, q, r are non-negative integers. We will focus on a particular class of Hirota bilinear 
equations in 2 + 1 dimensions: 

F{Dx, Dy, Dt)f ■ f = {DtP{Dx) + DyQ{Dx) + R{Dx))f ■ / = 0, (2.5) 

where P and Q are nonzero odd polynomials and i? is a nonzero even polynomial, namely, 
P, Q and R are nonzero polynomials and satisfy 

P{-z) = -P{z), Qi-z) = -Qiz), R{-z) = R{z). (2.6) 

When the Hirota operators act on exponential functions, the following derivative for- 
mula holds: 

DlDlDle"^^ ■ e'^^ = {h - k2f{li - hYiooi - ^2)"e''^+''^ (2.7) 



where ijj = kjX + ljy + Ldjt + rjjo, j = 1,2, with kj, lj,ujj, tjjq being constants. More generally, 
we have 

G{D,, Dy, Dt)e^^ ■ e^' = G{k, - k^, h - h, uji - a;2)e'^^+''% (2.8) 

where G is a polynomial in the three variables. This derivative formula will be a crucial 
key to our success in generating one-periodic and two-periodic wave solutions. 

We would like to consider the multi-dimensional special Riemann theta function solu- 
tion [28]: 

/ = f(x, y,t)= J2 e2"*<'''">+"*<™''^\ (2.9) 

where (■, ■) is the standard inner product of M^, n = (ni, ■ ■ ■ , un)'^ , rj = (r^i, ■ ■ ■ , 1]^)'^ 
with rjj = kjX + Ijy + Ujt + r]jQ, and r = (rpg)7VxAf is a symmetric matrix whose imaginary 
part is positive definite (i.e., Imr > 0). Base on (12.81) . we can compute G{Dx, Dt, ■ ■ ■) f ■ f 
generally for such a Riemann theta function / [29] , but we will make direct computations 
to provide a complete solution process and capture more of special solution structures. 

2.1 One-periodic wave solutions 

Let us first consider the case of = 1. Then the Riemann theta function in (12. 9p becomes 

oo 

f = f{x,y,t)= J2 e2™''+™'^ (2.10) 

n=— oo 

where Imr > and rj = kx + ly + ut + rjo with k, I, lo, rjo being real constants. 
Based on the derivative formula (12. 8p . we can compute that 



oo oo 
.2, 



n=—OD m=—oo 

oo oo 

n=—oo m=—oo 
oo oo 

= ^ ^ F(2^i(n-m)A:,27ri(n-m)/,27ri(n-m)o;)e2"*("+")^+"^("'+"')" 

n=— oo m=— oo 
oo oo 

= { ^ F(2^2(2n-m')fc,27ri(2n-m')/,27rz(2n-m')a;)e"[("'+(''-™')']"} 

%'=—oo n=~c 
oo 

F{m' 



m'=—oo n=~OD 
oo 



where the new summation m' = m + n has been introduced and F{m') is defined by 

oo 

F{m')= J2 i^(27r«(2n-m')A;,27ri(2n-m')/,27r?(2n-m»e"'["'+("-"')']". (2.11) 

n=— oo 

Shifting index n by introducing n' = r?, — 1, we have 

oo 

F(m') = i^(2^?(2n-m')fc,27rz(2n-m')/,27ri(2n-m')cj)e"^l"'+("-"^')'l" 



n=— oo 
oo 



F{2TTi[2n - (m' - 2 )]A;, 2Tci[2n' — {m! — 2)]/, 2Tii\2n' — {m! — 2)]a;) 



n — — oo 



^g7ri{n'2 + [n'-(m'-2)]2}Tg27ri(m'-l)7 



= F{m' - 2)e2-^(™'-i)-, m' G Z. 
It then follows that if F(0) = F(l) = 0, then F(m') = for all m' G Z. 

Noticing the specific form of the equation (12.51) . one-periodic wave solutions can be 
obtained, if we require 

oo 

F(0) = Y ['inniujP{4mrik) + 4mvilQ{4mrik) + R{Anmk)y'^^^" = 0, 



n=— oo 
oo 



F(l) = Y [2(2^ - l)mujP{2{2n - l)7rik) + 2{2n - \)i:HQ{2{2n - \)'Kik) 

n=— oo 

+i?(2(2n - l)7riA;)]e(2"'-2"+iK^^ = o. 



(2.12) 



Upon introducing 



an (A;) = ^ \n'KiP{4n'Kik)e^'^''^'\ 

n=—oo 
oo 

a^^{k) = J2 4n^iQ(4n^^A;)e2"'"*", 

n=— oo 
oo 

a2i{k)= Y 2(2n- l)7rzP(2(2n- l)7riA:)e(2'''-2''+^)"^", 

n=— oo 
oo 

a22{k) = J2 2(2Ti - l)TTiQ{2{2n - l)7r«fc)e(^"'-^"+^)"*", 



(2.13) 



and 



b^{k) = - Y Ri^nnik) 



n=— oo 

oo 



(2.14) 



b2{k) = - Y R{'2{2n - l)nik)t 



,(2n2-2n+l)7rir 



the linear system (12.121) of uj and / can be compactly written as 



aii(A;)u; + ai2(/i;)/ = 6i(A;), a2i{k)uj + a22{k)l = h2{k) . (2.15) 

We will see that there are a lot of choices for the angular wave number k. In order to 
generate real solutions (cu,/) to the system f l2.15p . we assume that 

Rer = 0. (2.16) 

The determinant of the coefficient matrix A{k) = {ars{k))2x2 is a polynomial in k, and so, 
if det{A{k)) ^ (this condition will be satisfied in our concrete examples), then 

Ao:= {k em\det{A{k)) = 0} (2.17) 

is either an empty set or a finite set. This guarantees the existence of real solutions {u, I) to 
the system f l2.15p at least for k ^ Aq. About nonzero solutions, we can have the following 
analysis. 

If deg(i?) = 0, i.e., R = c, where c is a nonzero real constant, then it follows from 
f l2.14p that b{k) does not depend on k and 

bik) = {hik),b2{k)f ^0, 

and so, there is the unique nonzero solution of {uj, I) to the system (12.150 for k ^ Aq. 
If deg(i?) > 2, then 

Bo := {k e M |(6i(A:))2 + b2ik)Y = 0} (2.18) 

is either an empty set or a finite set, since each of bi{k) and 62 (A;) is a polynomial in k 
of degree deg(i?). Therefore, there is the unique nonzero solution of {ujJ) to the system 
for k^AoU Bo. 
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2.2 Two-periodic wave solutions 

Let us second consider the case of = 2 and the corresponding two-periodic wave solutions. 
Similarly, based on the derivative formula fl2.8p and introducing m' = n + m, we can have 

F{D^, Dy, Dt)f ■ f = F{D^, Dy, /}^)e2^^('?'")+^^(™'") . e2^iM+ni{rm,m.) 

= F{27ri{n - m, k),2m{n - m, l),2m{n - ^))e2-^(^'"+™)+-*((™."^)+(™'")) 

= Y -^(27r«(2r2-m',A;),27r2(2n-m',/),27r«(2n-m',u;)) 

^ ■!Ti((T(n—m'),n—Tn') + (Tn,n))2iTi(r],m') 

Ac C 



2iTi{r],m') 
2)^ ' 



where F{m[,m'2) = F{m') is defined by 

F(m;,m^) = Y i^(27^^(2n-m^A:),27r^(2r^-m^/),27r^(2n-m^w))e"*(<"("-'"')"'>+<™'^^ 

(2.19) 

Shifting index n as n' = n — er with r = 1 or r = 2, where Ci = (1, 0)"^ and 62 = (0, 1)"^, 
we can compute that 

F{m[,m'2 - 2)e^''*('"2"^)""22+2Wiri2^ ^ ^ 2, 
where r = {rpq) 2x2- It now follows that if 

F(0, 0) = F(0, 1) = 0) = 1) = 0, (2.20) 
then F{m[, 1712) = for all m[,m2 G Z. 

For our selected equation (12. Sp . we have 

F(mi, 7712) = [27rz(2n — m, u)P{2'iri{2n — m, k)) 



n£2 



+2ni{2n - m, l)Q{2iTi{2n - m, k)) 
+R{2in{2n - m, A;))]e"^«"("-'")">+<™''^», 

where we set 



(2.21) 



m= (mi, 771,2)^, n=(ni,?7,2)^, k = {ki,k2)^, ^=(^1,^2)^, u = {lji,lj2)^ ■ (2.22) 



For simplicity, define 

where m*^'"'' = {m^[\ mg^'*)^, 1 < r < 4, are given by 

m« = (0,0f , m(2) = (0,lf , m(3) = (l,of , m^^) = (1, if . 
Then, upon introducing 

oo 

ars{k) = 2ni{2ns-mP)P{2ni{2n-m^''\k))er{n), 

ni ,712=— oo 
oo 

ar,s+2{k) = 27r2(2r2, - m'^'''>)Q{27ri{2n - m^^^, 

ni ,n2=— oo 

where 1 < r < 4 and 1 < s < 2, and 

oo 

b^[k) = - Y R{2TTi{2n-m''''\k))9r{n), 1 <r <4:, 

ni,n,2=— oo 

the hnear system (12.201) of (a;, I) can be compactly written as 



(2.23) 



(2.24) 



(2.25) 







- h{k) - 




= b{k) = 


b2{k) 


h 


h{k) 






. Hk) _ 



(2.26) 



where A{k) = (ars(^))4x4- If 't is purely imaginary, i.e., it satisfies (12.161) . then A{k) and 
b{k) are real, due to our assumption on the polynomials P, Q, R. Note that if det(y4.(A;)) ^ 
(this condition will be satisfied in our concrete examples), then 

Ao:= {k eM.^\det{A{k)) =0} (2.27) 

is either an empty set or a finite set. 

Now if deg(i?) = 0, i.e., R = c, where c is a nonzero real constant, then it follows from 
(I2.25P that b{k) does not depend on k and b{k) ^ 0, and so, there is the unique nonzero 
solution of {001,002, h, h) to the system fl2.26p for k ^ Aq. 



If deg(i?) > 2, then 



Bo:= {keR^\J2ibr{k)y = 0} 



(2.28) 



r=0 



is either an empty set or a finite set, since each of br{k), 1 < r < 4, is a polynomial in ki 
and k2 of degree deg R. Therefore, there is the unique nonzero solution of {lui,lu2, h, h) to 
the system ([226]) for k ^ AqU Bq. 
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3 Two illustrative examples 

Let us illustrate our idea of generating one-periodic and two-periodic wave solutions through 
two particular Hirota bilinear equations. The first example is 

Ut + U^Xy - 3UUy - 3Ua;f = 0, Vr, = Uy, (3.l) 

in the physical field. This nonlinear equation is related to the breaking soliton equation 

m- 

Ut + Uxxy - 4:UUy - 2Ur,d~'^Uy = 0, 

and it can be transformed into 

{D,D, + DyDl + c)f-f = 0, (3.2) 
where c can be an arbitrary function of y and t, under the transformation 

u = -2{\nfU, v = -2{\nf)^y. (3.3) 

Actually, we have 

Ut + U^xy - ^UUy - 3U^V = - [ — j^. 

The second example is 

Ut + Uxxxxy - {5uxxV + WuxyU - 15u'^v)x = 0, Vx = Uy, (3.4) 

in the physical field. This nonlinear equation can be transformed into 

{D,D, + DyDl + c)f-f = 0, (3.5) 

where c can be an arbitrary function of y and t, under the same transformation (13.31) . 
Similarly, we have 



The involved arbitrary function c of y and t shows the diversity of solutions to (2+1)- 
dimensional differential equations. 

To generate one-periodic and two-periodic wave solutions by the solution method in the 
last section, we need to assume that the above function c is constant, based on which the 
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angular wave number I (or the pair of angular wave numbers li and I2) and the frequency 
00 (or the pair of frequencies ui and 0^2) are constant and thus the derivative formula f l2.8p 
will hold. Obviously, we have 

P{z) = z, Q{z) = z^ R{z) = c, (3.6) 

for the equation (13.11) and 

P{z) = z, Q{z) = z\ R{z) = c, (3.7) 

for the equation (13. 4p . The polynomials P and Q defined above are odd and the polynomials 
R defined above are even, and so, the property (12. 6p is satisfied. The determinants of the 
corresponding coefficient matrices of the linear systems (12.151) and (I2.26p are not identically 
equal to zero, namely, 

det{A{k)) ^ and det(A(A;i, A^a)) ^ 
in the two examples. For instance, in the case of one-periodic wave solutions, we have 

det{A{k)) = ak^ or det{A{k)) = bk^, (3.8) 

where 

00 00 
a = -2567r6 Yl ^'e^-'-- ^ (2n - i)4e(2n^-2-+i)-- 

71=— 00 n=—oo 
00 00 

+10247r6 n^e^"'-- ^ (2n - l)2e(2"'-2"+i)--, 

n=— 00 n=—oo 
00 00 

b = 10247r« Yl ^'e^"'"'" Y ' l)^e(2"'-2"+i)-- 

n=— 00 n=— 00 

00 00 

-163847r8 Y "'e^"'-- Y (2^ ' l)2e(2"'-2-+i)--. 

n=— 00 n=— 00 

A direct computation by Maple 11 with Digits = 30 shows that 

a\r=o.u ~ 4563.212514, a|^=o.2i ~ 140396.7042, a|^=o.5i ~ 25831.08621, 
b\r=o.ii ~ 11012599.24, b\r=o.2i ~ 28544399.95, 6|^=o.5i ~ -4884657.870, 

which are all nonzero. Generally, our general analysis made before is valid for the two 
equations (13.21) and (13. 5p . and so, one-periodic and two-periodic wave solutions to the two 
(2 + l)-dimensional nonlinear equations possessing Hirota bilinear forms, (13.10 and (13.40 . 
can be computed explicitly. 
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4 Conclusion and remarks 



The Riemann theta functions have been used to generate one-periodic and two-periodic 
wave solutions of a particular class of (2+l)-dimensional Hirota bilinear equations, and the 
corresponding solution analysis has been made to guarantee the existence of such multi- 
periodic wave solutions. Two illustrative examples: 

Ut + Uxxy — 'iuuy — 3uxV = and Ut + Uxxxxy — i^UxxV + lOuxyU — 15u^v)x = 0, 

where Vx = Uy, have been discussed in details, along with their one-periodic and two- 
periodic wave solutions involving an arbitrary purely imaginary Riemann matrix. 

Our solution analysis provides a way to construct one-periodic and two-periodic wave 
solutions to (2+l)-dimensional nonlinear differential equations. It allows different angular 
wave numbers k (or different pairs of angular wave numbers ki and /C2), but the angular 
wave number / (or the pair of angular wave numbers li and I2) and the frequency u (or the 
pair of frequencies Ui and 002) are determined in terms of each angular wave number k (or 
each pair of angular wave numbers ki and k2) and hence the obtained solutions describe 
one-dimensional propagation of waves. 

We also remark that the proposed approach can be applied to other nonlinear differ- 
ential equations. For example, the following combined equation with the Sawada-Kotera 
vector field: 

Ut + Uxxy - 3uUy - 3UxV + Uxxxxx - l5{uUxx " U^)x = 0, Vx = Uy, 

can be analyzed similarly. Under the transformation (13.31) . this equation can be put into 
the following bilinear equation: 

{DtDx + DyDl + Dl + c)f-f = 0, 

where c can be an arbitrary function of y and t. The corresponding polynomials P,Q,R 
read 

P{z) = z, Q{z) = z\ R{z) = z^ + c, (4.1) 

where c is assumed to be constant. Therefore, the same analysis on one-periodic and 
two-periodic wave solutions will work for this equation as well. On the other hand, soliton 
solutions to the equations (13. ip and (13.41) can be computed by using Hirota's direct method. 
For example, one soliton solutions to the equations (13.11) and (13. 4p are determined by 

f = I ^ e^kH+kx^ky J ^ I ^±kH+kxTky ^ A: - arbitrary coust., 
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respectively. This can also be verified by using (12.81) . It should be, however, interesting to 
establish any relations between soliton solutions and multi-periodic wave solutions. 

It is our hope that our analysis on one-periodic and two-periodic wave solutions made 
for the particularly selected class of Hirota bilinear equations could help to better under- 
stand the diversity and integrability of nonlinear differential equations. 
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